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Abstract. We calculate three jet cross sections in photoproduction using exact matrix elements for the 
direct and resolved contributions. Numerical distributions are presented in a generic, irreducible set of 
variables that allows to disentangle the dynamics of partonic QCD subprocesses from each other and from 
pure phase space distributions. The results are compared to preliminary data from the ZEUS collaboration 
at HERA. It is found that the largest contribution comes from photon-gluon fusion in the mass range 
36 GeV < Ms-jet < 80 GeV. The measured leading jet scattering angle distribution is consistent with 
the (-channel exchange of a massless fermion in 2 —> 2 scattering, where the third parton is assumed to 
arise from soft bremsstrahlung. The data are inconsistent with pure phase space and Rutherford scattering 
distributions. 


1 Introduction 

Photoproduction of jets has been studied to a great ex¬ 
tent at HERA since the DESY electron-proton collider 
turned on in 1992 by the HI |l[ and ZEUS J 2 J collabora¬ 
tions. Until recently, however, the limited luminosity only 
allowed for measuring inclusive single jet and dijet cross 
sections. In general, the data have shown good agreement 
with QCD predictions accurate to next-to-leading order 
in perturbation theory || . They have confirmed the exis¬ 
tence of direct and resolved partonic contributions for real 
and slightly off-shell photons and begin to show sensitiv¬ 
ity to different parametrizations of the parton densities in 
the photon provided that uncertainties from the underly- 
inff hadronic event and the jet definition are under control 

li- 

On the other hand, multijet production has been mea¬ 
sured some time ago in purely hadronic collisions at the 
CERN SPS collider by the UA1 collaborationfl and at 
the Fermilab Tevatron collider by the CDF |[7J and DO 
H collaborations. Its importance lies not only in testing 
perturbative QCD to higher orders, but also in the search 
for new phenomena. Many analyses for the production of 
Standard Model and new particles, e.g. the top quark, the 
Higgs boson, or supersymmetric particles, rely on many 
hard jets in the final state. QCD multijet production then 
is a significant background for these searches and has to 
be well understood. 

At HERA the single jet transverse energy and dijet 
mass distributions in photoproduction have shown no de¬ 
viations from the next-to-leading order QCD predictions 
so far. However, higher integrated luminosity has now per¬ 
mitted ZEUS to produce the first three jet analysis Q. 
Even a limited number of four jet events has been ob¬ 
served. It is therefore interesting to compare these new 
data with QCD predictions and look for deviations as sig¬ 


nals of new physics. Since next-to-leading order cross sec¬ 
tions for the photoproduction of three jets are not avail¬ 
able, such an analysis is still restricted to leading order in 
perturbation theory and will therefore be subject to large 
renormalization and factorization scale uncertainties. As 
a possible way out, we will concentrate in this paper on 
normalized distributions that are largely independent of 
both higher order corrections and scale choices as are the 
shapes of the distributions. We briefly review the leading 
order cross section formalism and the relevant kinematic 
variables before we present numerical results for both dy¬ 
namic QCD and pure kinematic phase space distributions. 
We disentangle the different partonic subprocesses that 
contribute to photoproduction of three jets and show that 
the leading jet scattering angle distribution in three jet 
production is closely related to two jet production. 

The paper is organized as follows: In Sect. |^, we briefly 
review the leading order 2 —> 3 parton scattering pro¬ 
cesses. In Sect. H we define the relevant phase space vari¬ 
ables and the hadronic three jet cross sections. Our numer¬ 
ical results are presented in Sect. |], and we summarize our 
analysis in Sect. ||. 


2 Leading order 2 —> 3 parton scattering 

Our predictions for three jet cross sections in photopro¬ 
duction are based on the leading order 2 —> 3 matrix el¬ 
ements |AT | 2 averaged over initial and summed over final 
spin and color states. The partonic cross section can then 
be calculated through 

°i f = ^l- A/ f| 2 ( 27r ) 4 ^(Pi+P2-P3-P4-P5)n^I <5 ^i)- 

(1) 
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Fig. 1. Leading order 2^3 Feynman diagrams for direct 
photoproduction. 


We denote the four-momenta of the incoming and outgo¬ 
ing partons by £> 1,2 and £> 3 , 4,5 and the partonic center-of- 
mass energy squared by s = (p\ + P 2 ) 2 ■ 

Both direct and resolved processes have to be taken 
into account. The generic diagrams for direct 2 —» 3 scat¬ 
tering, where the photon interacts directly with a parton 
in the proton, are displayed in Fig. These diagrams arise 
from the leading order 2 —> 2 QCD Compton scattering 
process 


\M\\ 


16Tr 2 aa s 2Cp 



( 2 ) 


through the emission of an additional gluon (left diagram 
in Fig. [lj) or the splitting of the final state gluon into a 
quark-antiquark pair (right diagram in Fig. |lj). a and a s 
denote the electromagnetic and strong coupling constants, 
s, t, and u are the usual Mandelstam variables for 2 —> 2 
scattering, and Cf = (Nq — l)/(2 Nc) with Nc = 3 is the 
SU(3) color factor. It is worth noting that both the QCD 
Compton and the crossed photon-gluon fusion processes 
proceed through a massless fermion exchange in the t- 
cliannel, which leads to a single 1/t pole in the matrix 
element. Both diagrams in Fig. [I] lead to an extra factor 
of a s , when the matrix elements are squared, so that the 
direct 2 —>■ 3 cross sections are of G(aa 2 ). The outgoing 
quark-antiquark pair can be of the same or different flavor 
than the incoming quark. The diagrams for incoming anti¬ 
quarks or gluons can be obtained from Fig. [j] by crossing 
a final quark or gluon line with the incoming quark line. 
All possible topologies and orders of outgoing particles 
have to be considered for the full matrix elements of the 
processes, although they are not shown here explicitly. The 
complete result in d dimensions can be found in [ 10 ] in a 
very compact notation. Since there are no singularities 
present in three jet production, we can simply set d = 4. 

If the photon resolves into its partonic content before 
the hard scattering takes place, the effective hard scat¬ 
tering is of purely partonic nature. The relevant generic 
diagrams are displayed in Fig. I They are of 0(a ]!) and 
arise from the underlying 2 —> 2 parton-parton scattering 
processes through radiation of an additional gluon in the 
final state. Most of the 2 —> 2 parton-parton processes 
proceed through the exchange of a massless vector boson 
in the t-channel leading to a double pole 1 /t 2 in the matrix 
elements. For example, the matrix element for the process 
qq' —> qq' is given by 


m 


qq'^qq' 


= 16n 2 a 


2 

s 


Cf s 2 + u 2 

ivc t 2 


( 3 ) 
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Fig. 2. Leading order 2 —> 3 Feynman diagrams for resolved 
photoproduction. 


Crossing one or two outgoing gluon lines into the initial 
state in Fig. leads to the gluon initiated processes with 
quarks in the final state, and crossing an incoming and 
outgoing quark line leads to the processes with incoming 
anti-quarks. Process c) is symmetric under the interchange 
of the three gluons and under the interchange of the quark 
and anti-quark. Process d) is completely symmetric under 
the interchange of any of the five gluons. The complete 
result for the 2 —> 3 parton scattering matrix elements in 
d dimensions can be found in fll] ] , where we can set d = 4 
again. 


3 Hadronic three jet cross sections 

At HERA, positrons and protons are currently collided 
with energies of E e = 27.5 GeV and E p = 820 GeV. The 
hadronic three jet cross section can be calculated from the 
partonic cross section off 5 with the help of the factoriza¬ 
tion theorem, 

a ep 3 et ( s H) = J dx e f e (x e ,M 2 )dx p fp(Xp,M 2 ) 

alf{x e x p SH,M 2 ). (4) 

M is the factorization scale and y / J]j = 1 /( p e + Pp) 2 = 
300 GeV is the hadronic center-of-mass energy. f e denotes 
the parton density in the positron and is given by a con¬ 
volution of the density of photons in the positron with the 
density of partons in the photon, 

f e (x e ,M 2 )= J ^f:(y)f^(x e /y,M 2 ). (5) 

fjiy) is given by the Weizsacker-Williams-Approximation 
[12], and we include the subleading terms calculated in 
[13|. Experimentally the photon has a virtuality of Q < 1 
GeV and carries a fraction 0.2 < y < 0.8 of the positron 
energy. For the parton densities in the photon f^' q ’ 9 and 
proton/®’ 9 , we use the leading order parametrizations by 
GRV JL4| and CTEQ4 ]li| except where stated otherwise. 
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Asymptotically, the parton densities in the photon are of 
0 {a/a s ), so when they are convoluted with the resolved 
2 —> 3 matrix elements of 0(o? s ), the hadronic cross sec¬ 
tion is of the same order as the direct contribution. The 
strong coupling constant a s (/z) is calculated in the one- 
loop approximation with five flavors and = 181 MeV 
as found in the global fit of the parton densities in the 
proton. The renormalization and factorization scales are 
set to the maximum transverse jet energy in the event, 
// = M = inaxf /iT i, At.2- A'7- 3). 

In the ZEUS three jet analysis, jets are defined by the 
KTCLUS algorithm in inclusive mode |l(|. This defini¬ 
tion uses a parameter R = 1 and combines two hadronic 
clusters i and j if 

dij = min (E^, E$ tj )[(rn - r]j ) 2 + (fc - (/>j) 2 ]/R 2 (6) 

is smaller than 

di=E^. (7) 

The transverse energy Et, pseudorapidity 77 , and azimuthal 
angle <f> of the combined cluster are calculated from the 
transverse energy weighted sums of the two pre-clusters, 

Et = Er,i + Etj, (8) 

V = [ET,iVi + E T,jr]j}/ E T , (9) 

4> = [ ET,i<t>i + ET,i<j>j\/E t- ( 10 ) 

The longitudinal momentum fractions x p and x e of a par- 
ton in the proton and electron can be determined from 
the observed three jet final state through four momentum 
conservation 


5 

x P ,e = ^2 Erte^/y/sZ, 
i =3 


(ii) 


where the plus sign applies to the proton going into the 
positive z direction and the jet pseudorapidities rji are 
defined in the center-of-mass frame. The longitudinal mo¬ 
mentum fraction of a parton in the photon is given by 
x 1 = x e /y. The three jets are required to have transverse 
energies of Et, 3 ,Et ,4 > 6 GeV and Et ,5 > 5 GeV, and 
their pseudorapidities must lie within 1773 , 4 , 5 ! < 2.4 in the 
laboratory frame |fj. The KTCLUS algorithm is favored 
experimentally and in higher order theoretical predictions 
since it contains no overlap or double counting ambigu¬ 
ities. Leading order QCD predictions, however, lack the 
possibility to implement or distinguish a particular jet def¬ 
inition, since each parton is identified with a jet. 

For N massless jets, one can choose 3-/V —4 parameters 
that should span the multijet parameter space, facilitate 
a simple interpretation within QCD, and allow for a com¬ 
parison of the N — 1-jet to the IV-jet cross section. In 
the case of N = 3, the conventional choices are the three 
jet mass M 3 _j et , required by ZEUS to be larger than 50 
GeV, and four dimensionless parameters. The Dalitz en¬ 
ergy fractions 


Xi = 


2 Ej 

M 3 —jet 


( 12 ) 


specify how the available energy is shared between the 
three jets. They are ordered such that X3 > X4 > X5. Since 
X 3 + X 4 + X 5 = 2, only X 3 and X 4 are linearly independent. 
The third and fourth parameters are the cosine of the 
scattering angle between the leading jet and the average 
beam direction pav = Pi — P 2 , where the incoming parton 
1 is the one with the highest energy in the laboratory 
frame, 


cos 83 


PAVP 3 

IpavIIpsI ’ 


(13) 


and the angle between the three jet plane and the plane 
containing the leading jet and the beam direction, 


, Ps x pav P 4 x p 5 

c°s ^3 = 1 -n-r- 

P 3 X PAV P 4 X Ps 


(14) 


Instead of cos ^3 we will investigate distributions in ip 3 
itself. In the soft limit, where E 5 —> 0 and £ 3,4 —> 1, 
cos 83 approaches the 2 —> 2 center-of-mass scattering an¬ 
gle cos 8 * thus relating three jet to two jet cross sections, 
cos 8 * can be determined from the pseudorapidities of the 
two jets by 

cos 8 * = tanh |^ 7?3 ~ T?4 ^ (15) 

and is related to the 2 —> 2 Mandelstam variables by 
t = —is(l — cos 0 *) and u = — ^ s (l + cos 9*). (16) 


Of course, the third jet must not be too soft (or the hard 
jets not too hard) to avoid soft singularities that would 
have to be absorbed into the next-to-leading order dijet 
cross section. This is achieved by a cut on X 3 < 0.95. 
However, since the energy of a jet is always larger than 
its transverse energy, the cut on Et ,5 > 5 GeV already 
insures the absence of soft singularities. The three jets 
also have to be well separated in phase space from each 
other and from the incident beams to avoid initial and 
final state collinear singularities. This is insured by a cut 
on | cos 83 1 < 0.8, by the cuts in the pseudorapidities, and 
by the isolation cuts in the jet definition. 


4 Numerical results 

Having explained the relevant theoretical details in the 
last two sections, we can now turn to the presentation 
of numerical results for the photoproduction of three jets 
at HERA. We first look for discrepancies between the¬ 
ory and ZEUS data Q in the three jet mass distribu¬ 
tion in Fig. |j. We find that the total theoretical predic¬ 
tion (full curve) describes the data well in shape and nor¬ 
malization. The agreement in normalization is, however, 
to some degree coincidental since there is still an uncer¬ 
tainty of a factor of two from the variation of the scales 
/j, = M £ [0.5; 2.0] x max(U , 7 ', 3 , Et, 4 , Et#) (shaded band). 
This theoretical uncertainty is much bigger than the sta¬ 
tistical experimental error of about 5% at small M 3 _j e t 
indicating the need to implement next-to-leading order 
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Fig. 3. Total cross section (full curve) for the photoproduction of three jets as a function of the three jet mass M 3 _j e t. We also 
show the variation of the absolute normalization due to the uncertainty in the scale choice (shaded band) and the contributions 
from direct photons (dashed), quarks (dotted), and gluons (dot-dashed) in the photon. The ZEUS data || agree well with the 
QCD prediction in shape and normalization within the statistical error. 


corrections in the three jet cross section. Since these are 
not available yet, we will subsequently resort to normal¬ 
ized distributions that are largely independent of both 
higher order corrections and scale choices. In Fig. |] we 
also show the individual contributions to the total cross 
section by direct photons (dashed curve, contributes more 
than 50%), by quarks (dotted curve, contributes less than 
40%), and by gluons (dot-dashed curve, contributes less 
than 10%) in the photon. From the inset smaller figure it 
becomes clear that the quark-initiated process gives the 


largest contribution for M 3 _j et < 36 GeV. We also com¬ 
puted but do not show here the individual contributions 
from quarks and gluons in the proton and find the gluon 
to give the largest contribution for M 3 _j et < 80 GeV. In 
the region of 36 GeV < M 3 _j et < 80 GeV, the photon- 
gluon fusion process accounts for one third of the total 
cross section. 

In Fig. ||, we show the total three jet cross section 
per bin as a function of the longitudinal momentum frac¬ 
tions x 1 of the partons in the photon (left) and x p of 
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Fig. 4. Total three jet cross sections per bin (full curves) as functions of the longitudinal momentum fraction * 7 of a direct 
photon (dotted), quark (dot-dashed), or gluon (long-dashed) in the photon (left) and of the longitudinal momentum fraction 
x p of a quark (dotted) or gluon (dot-dashed) in the proton (right). Also shown are the pure phase space distributions (dashed 
curves) which differ clearly from the dynamical QCD distributions. 


the partons in the proton (right). These distributions re¬ 
flect the situation near the cut at M 3 _j et > 50 GeV, since 
the three jet cross section falls off exponentially with the 
three jet mass. In the left figure, the total cross section 
(full curve) is presented together with the direct peak at 
:r 7 = 1 (dotted), which has been divided by the bin width, 
and the quark (dot-dashed) and gluon (long-dashed) ini¬ 
tiated processes. The kinematic cuts restrict the phase 
space to Xj > 0.1 (dashed). The quark initial state is big¬ 
ger than the gluon above a ; 7 = 0.3, and both contributions 
are larger (smaller) in the GRV parametrization than in 
GS96 (nj when x 1 > 0.35 (x 7 < 0.35). For the longitu¬ 
dinal momentum fraction x p of the parton in the proton 
(right side of Fig. ||), the kinematic cuts restrict the phase 
space to x p > 2.5 • 10 -2 and render contributions above 
x p > 0.5 insignificant. The gluon (quark) in the proton 
dominates for x p < 0.16 (x p > 0.16), and at low three 
jet masses, the gluon contributes about two thirds of the 
total cross section. 

We now turn to normalized distributions in the four 
dimensionless three jet parameters X3, X4, cos 6 * 3 , and 1/J3. 
Some of the theoretical predictions have already been in¬ 
cluded in the presentation of the experimental data [j9j. 
Fig. I shows the three jet cross section as a function of the 
energy fractions of the leading (left) and next-to-leading 
(right) jet x 3 and X4, normalized to the total cross section. 
The prediction from the O(aa^) QCD matrix elements 
(full curves) are rather similar to the pure phase space 
distributions (dashed curves), but the data slightly prefer 
the QCD predictions. The statistical error bars are still 
too large for a definite conclusion, however. We have also 
studied the relative contributions from the different par- 


tonic sub-channels. They are flat in a; 3 , but photon-quark 
scattering clearly dominates when X 4 ~ X 3 —> 1 . 

The distribution in the cosine of the scattering angle 
of the fastest jet cos 0 3 is displayed in the lower plot of 
Fig. H together with the distribution in the dijet scatter¬ 
ing angle cos 8 * (upper plot) measured earlier by ZEUS J|] 
and compared to theoretical predictions in ||. They are 
both of O(aa^) which is the leading order for the three 
jet case but next-to-leading order for the dijet case. The 
wiggles in the latter stem from limited statistics in the 
Monte Carlo integration. Of course, the cuts in the dijet 
and three jet cross sections differ slightly, but the most 
important cuts on the jet mass M 2 _j e t > 47 GeV and the 
scattering angle |cos 0 *| < 0.8 are in fact very similar. 
The differences will also mainly affect the normalization 
and should be insignificant in distributions normalized at 
cos 6 = 0 as they are presented here. Let us first con¬ 
centrate on the dijet case. As expected it is symmetrical 
in cos 9*. The data show good agreement with the QCD 
prediction (full curve) but clearly disagree with the pure 
phase space distribution (dashed curve). It might be sur¬ 
prising at first sight that the data also disagree with the 
Rutherford scattering form at small angle, (1— | cosd*|) -2 . 
On the other hand, the data agree very well with the less 
singular form (1 — | cos0*|) . This can be understood if 

we remind ourselves that Rutherford scattering is charac¬ 
teristic of the exchange of a massless vector boson in the 
f-channel 


\M\ 2 oct ~ 2 


— s(l — cos 0 *) 
2 v ' 


-2 


(17) 
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Fig. 5. Dependence of the total cross section (full curves) and phase space (dashed curves) on the energy fractions £3 (left) 
and X 4 (right) of the leading and next-to-leading jets. The ZEUS data [Iff show a slight preference of the full QCD prediction 
over phase space. 


as present in most of the resolved (parton-parton) scat¬ 
tering processes whereas both direct processes proceed 
through a massless fermion exchange in the f-channel with 
less singular behavior, 


dictions when three jet predictions were still unavailable 

m- 

To study deviations from the (1 — cos 9) 2 behavior it 


is convenient to plot the data in terms of the variable 119 


\M\ 2 oc t~ l 


- —s(l — cos 9*) 


(18) 


In addition, the process 75 —> qg has an s-channel con¬ 
tribution without any singular behavior. Therefore, the 
cos 9* distribution provides evidence that direct processes 
and photon-gluon fusion in particular are the most impor¬ 
tant subprocesses in this kinematic region. Turning now 
to the three jet cross section, we observe that both the 
0(aa 2 ) QCD (full curves) and phase space distributions 
(dashed curves) are asymmetric in cos 9 3 which motivates 
why we kept the sign in the abscissa and did not present 
curves in the absolute value |cos 0 | as one usually does 
in the dijet case. Again the data clearly prefer QCD (full) 
over phase space (dashed) and the single pole in t (dotted) 
over the double pole (dot-dashed). In fact, the dependence 
of the three jet cross section on the scattering angle of the 
fastest jet is strikingly similar to the dijet distribution in 
the center-of-mass scattering angle. The deviation in the 
backward region of the three jet case cos 9 3 < 0 , when the 
pseudorapidity of the fastest jet 773 < 774,5 is tending to¬ 
wards the direction of the photon beam, can be traced to 
a large extent to the asymmetric phase space. Since three 
jet production proceeds dominantly through single gluon 
bremsstrahlung, we can relate this result to the underly¬ 
ing 2 —> 2 process with the interpretation given there. In 
fact it is interesting to note that the CDF collaboration 
compared their first three jet data directly to two jet pre- 


1 + cos 6 
1 — cos 9 


(19) 


which is related to the experimentally measured pseudo¬ 
rapidities by 

X = e Vl ~ V2 (20) 

and removes the Rutherford singularity such that 

da . . 

— ~ constant ( 21 ) 

d X 

if one ignores the angular dependence and scaling viola¬ 
tion in the coupling constant and the parton densities. 
In Fig. we present distributions in the angular variable 
X of the fastest of three jets, normalized at x = 1 (he. 
cos 9 = 0). The data show again a clear distinction from 
both pure phase space and Rutherford scattering at small 
angle (1— | cos0 |) -2 but agree very well with 0(aa 2 ) QCD 
prediction and also with the form (1 —| cos 0|) _1 . The QCD 
prediction includes, of course, scaling violation both in 
the strong coupling and in the parton densities. Absence 
of scaling violations would result in a constant prediction 
and would clearly disagree with the data. As before, the 
full QCD prediction is dominated by direct photon scat¬ 
tering. Therefore the long-dashed curve lies close to the 
full curve, whereas the processes initiated by the quark in 
the photon (long-dash-dotted curve) give a steeper shape 
as expected. The gluon initiated processes (wide-dotted 
curve) have an even different shape. The total resolved 
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Fig. 6. Cross sections for the photoproduction of two (upper plot) and three (lower plot) jets as functions of the center-of-mass 
and fastest jet scattering angles normalized at cos 9 = 0. The ZEUS dijet jHJ and three jet data jo| are well described by the 
QCD predictions (full curves), but not by the pure phase space distributions (dashed curves), and they favor the single 1/t pole 
for massless fermion exchange (dotted curves) over the double 1/f 2 pole for massless boson exchange (dot-dashed curves) in the 
f-channel. 


contribution (not shown), i.e. the sum of quark and gluon 
initiated contributions, lies between the quark initiated 
and total curves and accounts for the remaining difference 
of the direct initiated and total cross section. 

Finally we investigate the dependence of the three jet 
cross section normalized to the total cross section on the 
angle between the three jet plane and the plane contain¬ 
ing the leading jet and the average beam direction in 
Fig. |[ The full QCD curve again agrees well with the data 


as do the contributions from direct photons (long-dashed 
curve) and quarks in the photon (long-dash-dotted curve), 
whereas the pure phase space (dashed curve) has a com¬ 
pletely different shape. The gluon in the photon (wide- 
dotted curve) gives a slightly different shape. 
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x 3 

Fig. 7. Dependence of the total three jet cross section on the angular variable \3- The ZEUS data [Oj] agree with the QCD 
prediction (full curve), which is dominated by the direct contribution (long-dashed), and with the single pole in t (dotted). They 
exclude the pure phase space (dashed), quark (long-dash-dotted) and gluon (wide-dotted) initiated distributions as well as the 
double pole in t (dot-dashed). 



Fig. 8. Dependence of the total three jet cross section (full curve), the pure phase space (dashed curve), and the individual 
contributions from the partons in the photon (long-dashed, long-dash-dotted, and wide-dotted curves) on the angle between the 
three jet plane and the plane containing the leading jet and the average beam direction. The ZEUS data [[)] rule out the pure 
phase space and gluon initiated distributions. 
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5 Summary 

In this paper we have presented three jet cross sections 
in photoproduction using exact leading order matrix ele¬ 
ments. The normalization and shape of the QCD predic¬ 
tions show good agreement with ZEUS data, although the 
normalization is subject to rather large scale uncertain¬ 
ties. None of the distributions can be explained by the 
three jet phase space alone. Of the different direct and 
resolved contributions, direct and in particular photon- 
gluon fusion processes are the most important subpro¬ 
cesses. The fastest-jet scattering angle distribution looks 
very similar to the dijet center-of-mass scattering angle 
distribution indicating that three jet production proceeds 
mainly through single bremsstrahlung. The data are not 
consistent with Rutherford scattering for the exchange of 
a massless boson but with the less singular form for the 
exchange of a massless fermion in the f-channel. 
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